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1 . 


Introduction 


The global hydrological cycle is central to climate system interactions and the key to 
understanding their behavior. Rainfall and its associated precipitation processes are a key link 
in the hydrologic cycle. Fresh water provided by tropical rainfall and its variability can exert a 
large impact upon the structure of the upper ocean layer. In addition, approximately two-thirds 
of the global rain falls in the Tropics, while the associated latent heat release accounts for about 
three-fourths of the total heat energy for the Earth's atmosphere (Riehl and Simpson 1979). 
Precipitation from convective cloud systems comprises a large portion of tropical heating and 
rainfall. Furthermore, the vertical distribution of convective latent-heat releases modulates 
large-scale tropical circulations (e.g., the 30-60-day intraseasonal oscillation - see Sui and Lau 
1988), which, in turn, impacts midlatitude weather through teleconnection patterns such as 
those associated with El Nino. Shifts in these global circulations can result in prolonged 
periods of droughts and floods, thereby exerting a tremendous impact upon the biosphere and 
human habitation. And yet, monthly rainfall over the tropical oceans is still not known within a 
factor of two over large (5-degrees latitude by 5-degrees longitude) areas (Simpson et al. 1988, 
1996). Hence, the Tropical Rainfall Measuring Mission (TRMM), a joint U.S. /Japan space 
project, can provide a more accurate measurement of rainfall as well as estimate the four- 
dimensional structure of diabatic heating over the global tropics. The distributions of rainfall 
and inferred heating can be used to advance our understanding of the global energy and water 
cycle. In addition, this information can be used for global circulation and climate models for 
testing and improving their parameterizations. 

Cloud resolving (or cumulus ensemble) models (CRMs) are one of the most important 
tools used to establish quantitative relationships between diabatic heating and rainfall. This is 
because latent heating is dominated by phase changes between water vapor and smali, cloud- 
sized particles, which can not be directly detected using remote sensing techniques (though 
some passive microwave frequencies do respond to path-integrated cloud water). The CRMs, 
however, explicitly simulate the conversion of cloud condensate into raindrops and various 
forms of precipitation ice. It is these different forms of precipitation that are most readily 
detected from space, and which ultimately reach the surface in the form of rain in the Tropics. 
In addition, the highest science priority identified in the Global Change Research Program 
(GCRP) is the role of clouds in climate and hydrological systems, which have been identified 
as being the most problematic issues facing global change studies. For this reason, the 
GEWEX (Global Energy and Water Cycle Experiment) formed the GCSS (GEWEX Cloud 


System Study), specifically for the purpose of studying such problems. CRMs were chosen as 
the primary approach (GCSS Science Plan 1993; Moncrieff et al. 1997). 

The first pioneering one dimensional cloud model was developed by Dr. J. Simpson in 
the 1960s. A two-dimensional anelastic model that filtered out sound waves was developed by 
Drs. Y. Ogura and N. Phillips. The models were used to study cloud development under the 
influence of the surrounding environment. The ID cloud model was used extensively to study 
the cloud seeding problem. In the late 1970's, four three-dimensional cloud models were 
developed (Wilhelmson 1974; Miller and Pearce 1974; Sommeria 1976; Clark 1979; Klemp 
and Wilhelmson 1978; Cotton and Tripoli, 1978; and Schlesinger 1975, 1978). The effect of 
model designs (i.e., slab vs axis-symmetric, and 2D vs 3D) on cloud development and liquid 
water content were the major foci in 70's. Also, the dynamics of midlatitude supercells, that 
are usually associated with tornados, was another major focus in the 70's. After GATE, cloud 
ensemble modeling was developed to study the collective feedback of clouds on the large-scale 
tropical environment with the aim of improving cumulus parameterization in large-scale models 
(i.e., Tao 1978; Soong and Tao 1980; Tao and Soong 1986; Lipps and Helmer 1986; Krueger 
1988). The effect of ice processes on cloud formation and development, stratiform rain 
processes and their relation to convective cells, and the effect of wind shear on squall line 
development were the other major areas of interest for cloud resolving models in the 1980's. 
The impact of radiative processes on cloud development was also investigated in the late 80's. 
In the 1 990’s, cloud resolving models were used to study multi-scale interactions, cloud 
chemistry interaction, idealized climate variations, and surface processes. The cloud resolving 
model was also used for the development and improvement of satellite rainfall retrieval 
algorithms. Table 1 lists the major foci and some (not all) of the key contributors to cloud 
resolving model development over the past four decades. 

During the past 20 years, observational data on atmospheric convection has been 
accumulated from measurements by various means, including radars, instrumented aircrafts, 
satellites, and rawinsondes in special field observations (e.g., GATE, PRE-STORM, 
COHMEX, TAMEX, EMEX, TOGA COARE 1 and several others). This has made it possible 


1 GATE stands for GARP (Global Atmospheric Research Program) Atlantic Tropical 

Experiment, TAMEX for Taiwan Area Mesoscale Experiment, EMEX for Equatorial Mesoscale 
Experiment, PRE-STORM for Preliminary Regional Experiment for Storm Central, COHMEX for 
Cooperative Huntsville Meteorological Experiment, and TOGA COARE for Tropical Oceans Global 
Atmosphere (TOGA) - Coupled Ocean Atmosphere Response Experiment (COARE). 



for cloud resolving modelers to test their simulations against observations, and thereby 
improve their models. In turn, the models have provided a necessary framework for relating 
the fragmentary observations and helping to understand the complex physical processes 
interacting in atmospheric convective systems, for which observations alone still cannot 
provide a dynamically consistent four-dimensional picture. The past decades have also seen 
substantial advances in the numerical modeling of convective clouds and mesoscale convective 
systems (e.g., squall-type and non-squall-type convective systems), which have substantially 
elucidated complex dynamical cloud-environment interactions in the presence of varying 
vertical wind shear. With the advent of powerful scientific computers, many important and 
complex processes (which require extensive computations), such as ice-microphysics and 
radiative transfer, can now be simulated to a useful (but still oversimplified) degree in these 
numerical cloud models. Table 2 lists the key developments in the cloud resolving mode! 
approach for studying tropical convection over the past two decades. As shown, over the last 
20 years, these models have become increasingly sophisticated through the introduction of 
sophisticated (bulk-type) microphysical processes, radiation and boundary-layer effects, and 
improved turbulent parameterizations for subgrid-scale processes. In addition, an exponentially 
increasing computer resource has resulted in time integrations increasing from hours to days, 
domain grids boxes (points) increasing from less than 2000 to more than 2,500,000, and 3-D 
models becoming increasingly prevalent. The CRM is now at a stage where it can provide 
reasonably accurate statistical information of the sub-grid, cloud-resolving processes now 
poorly parameterized in climate models and 
numerical prediction models. 

2 . Goddard Cumulus Ensemble (GCE) Model 

The Goddard Cumulus Ensemble (GCE) model is a cloud resolving model, and its main 
features have been published by Tao and Simpson (1993) and Simpson and Tao (1993). The 
model is nonhydrostatic and model variables include horizontal and vertical velocities, potential 
temperature, perturbation pressure, turbulent kinetic energy, and mixing ratios of all water 
phases (vapor, liquid, and ice). The cloud microphysics includes a parameterized Kessler-type 
two-category liquid water scheme (cloud water and rain), and a three-category ice-phase 
scheme (cloud ice, snow and hail/graupel) mainly based on Lin et al. (1983) and Rutledge and 
Hobbs (1984). The Goddard microphysics scheme has several minor modifications, however. 
The first modification is the option to choose either graupel or hail as the third class of ice 
(McCumber et al. 1991). Graupel has a low density and a large intercept (i.e., high number 
concentration). In contrast, hail has a high density and a small intercept (i.e. low number 


concentration). These differences can affect not only the description of the hydrometeor 
population, but also the relative importance of the microphysical-dynamical-radiative 
processes. Second, a saturation technique was implemented by Tao et al. (1989b). This 
saturation technique is basically designed to ensure that supersaturation (subsaturation) cannot 
exist at a grid point that is clear (cloudy). This saturation technique is one of the last 
microphysical processes to be computed. It is only done prior to evaluating the evaporation of 
rain and snow/graupel/hail deposition of sublimation. A third difference is that all 
microphysical processes (transfer rates from one type of hydrometeor to another) are calculated 
based on one thermodynamic state . This ensures that all processes are treated equally. The 
opposite approach is to have one particular process calculated first modifying the temperature 
and water vapor content (i.e., through latent heat release) before the second process is 
computed. The fourth difference is that the sum of all the sink processes associated with one 
species will not exceed its mass. This ensures that the water budget will be balanced in the 
microphysical calculations. 

The following major improvements have been made to the model during the past seven 
year period: (i) The implementation of a multi-dimensional Positive Definite Advection 
Transport Algorithm (MPDATA, Smolarkiewicz and Grabowski 1990). All scalar variables 
(potential temperature, water vapor, turbulence coefficient and all hydrometeor classes) use 
forward time differencing and the MPDATA for advection. The dynamic variables, u, v and w, 
use a second-order accurate advection scheme and a leapfrog time integration (kinetic energy 
semi-conserving method), (ii) The development of an improved four-class, multiple-moment, 
multiple-phase ice scheme (Ferrier 1994), which resulted in improved agreement with observed 
radar and hydrometeor structures for convective systems simulated in different geographic 
locations without the need for adjusting coefficients (Ferrier et al. 1995). (iii) The inclusion 
of solar and infrared radiative transfer processes, which have been used to study the impact of 
radiation upon the development of clouds and precipitation (Tao et al. 1991, 1996) and upon 
the diurnal variation of rainfall (Tao et al. 1996; Sui et al. 1998) for tropical and midlatitude 
squall systems, (iv) The incorporation of land and ocean surface processes to investigate their 
impact upon the intensity and development of organized convective systems (Wang et al. 1996; 
Lynn et al. 1998). Mesoscale circulations, which formed in response to landscape 
heterogeneities represented by a land surface model, were crucial in the initiation and 
organization of the convection. 

A stretched vertical coordinate (height increments from 40 to 1150 m) is used to 
maximize resolution in the lowest levels of the model. Typically, a total of 1024 grid points are 



used in the horizontal with 500- 1000m resolution in the two-dimensional version of the GCE 
model. In the three-dimensional version of the GCE model, the horizontal resolution is usually 
2000m with 200 by 200 grid points. The time step is 5 to 10 s. Table 3 lists the characteristics 
of the GCE model. 

3 . Applications of the GCE Model to the Study of Precipitation Processes 

The application of the GCE model to the study of precipitation processes can be generalized 
into fourteen categories (Table 4). It has been used to provide essential insights into the 
interactions of clouds with each other (Tao and Simpson 1984, 1989a), with their 
surroundings, and their associated heat, moisture, momentum, mass and water budgets (Tao 
1978; Soong and Tao 1980, 1984; Tao and Soong 1986; Tao, Simpson and Soong 1987; Tao 
and Simpson 1989b), with radiative transfer processes (Tao et al. 1991, 1993a, 1996; Sui et 
al. 1998), with ocean surfaces (Tao et al. 1991; Wang et al. 1996, 2000), with idealized 
climate variations (Lau et al. 1993, 1994; Sui et al. 1994; Tao et al. 1999), and cloud draft 
structure and trace gas transport (Scala et al. 1990; Pickering et al. 1992; and a review by 
Thompson et al. 1997) and precipitation efficiency (Ferrier et al. 1996). The GCE model has 
also been used to convert the radiances received by cloud-observing microwave radiometers 
into predicted rainfall rates (Simpson et al. 1988, and a review by Simpson et al. 1996). 
Remote sensing of cloud-top properties by high-flying aircraft bearing microwave and other 
instruments is now beginning to provide powerful tests of the GCE model, particularly when 
such observations are augmented by simultaneous ground-based radar measurements (Adler et 
al. 1991 ; Prasad et al. 1995; Yeh et al. 1995). The GCE model has also been used to study the 
distribution of rainfall and inferred heating (Tao et al. 1990, 1993b, 2000a and b). In this 
paper, a brief discussion about the application of the GCE model to (1) cloud interaction and 
mergers, (2) convective and stratiform interaction, (3) mechanisms of cloud-radiation 
interaction, (4) latent heating profiles and TRMM, and (5) responses of deep cloud systems 
to large-scale processes will be provided. Comparisons between the GCE model’s results, 
other cloud resolving model results and observations will also be examined. 

3.1 Cloud Interactions and Mergers 

Field experiment data (e.g., FACE, Florida Area Cumulus Experiment; GATE, GARP Atlantic 
Tropical Experiment; and ITEX, Island Thunderstorm Experiment) has shown that the merging 
of shower clouds is a crucial factor in the development of organized convective complexes 
which are the major producers of rainfall in the tropics (Houze and Chang 1977), in the Florida 
peninsula (Simpson et al. 1980) and in the Maritime Continent region north of Darwin, 



Australia (Simpson et al. 1993). The observational data consisted of the calibrated radar and 
rain gauges. The mergers usually yield more than an order of magnitude more precipitation 
than unmerged cells. For example, Simpson et al . (1980) found that mergers were responsible 
for 86% of the rainfall observed, even though 90% of the cells were unmerged. Most of the 
increase in total rainfall comes from the increased areal extent and duration of the second-order 
mergers. [A first-order merger is identified as a consolidation of two or more previously 
independent single cell echoes, while a second-order merger is the result of the juncture of two 
or more first-order merged echoes (Westcott 1984).] 

However, the physical mechanisms which effect the merging process are not clearly 
specified through observational studies, largely because of the difficulty of measuring the air 
circulations in and around cumulus clouds. Westcott (1984) reviewed observational analyses 
of mergers in detail, and also raised some key questions concerning the mechanisms involved. 
From observational studies, several processes have been proposed as important in merging 
events. These processes fall into two main categories. The first involves addition of moisture 
to neighboring air, thereby reducing dilution by entrainment (Byers and Braham 1949; Scorer 
and Ludlam 1953; Malkus 1954). Moistening of the cloud environment can be accomplished in 
several ways. One source of moisture is precipitation falling from an overhanging canopy 
which produces a favorable environment for new convective growth. Dissipation of previous 
and nearby clouds also provides a moister, more favorable environment. The merging cells 
can be better protected from the entrainment of dry environmental air (Lopez 1978). The 
second category involves dynamic processes which enhance low-level convergence leading to 
new growth and merging. Low-level convergence can be enhanced by (1) collision of 
downdraft outflows (Simpson 1980; Simpson et al. 1980); (2) differential motions of cloud 
masses (Holle and Maier 1980; Cunning et al. 1982; LeMone 1989); and (3) hydrostatic and 
non-hydrostatic pressures response within the boundary layer (Cunning and Demaria 1986; 
LeMone et al. 1988). 

3.1.1 The GCE model simulation results 

A two-dimensional version of the GCE model was used with a GATE data set to study cloud 
interactions and merging (Tao and Simpson 1984). Over two hundred groups of cloud 
systems with a life history of over sixty minutes were generated under the influence of different 
combinations of the stratification and large-scale forcing (through a total of 48 numerical 
experiments). The GCE model results demonstrated the increase in convective activity and in 
amount of precipitation with increased intensity of large-scale forcing (lifting). In the GCE 



model simulation, a cloud merger is defined as a joining of the surface rainfall contour of 1 mm 
h'l . Additional criteria are also considered. The merged clouds need to join for at least 15 
minutes and the distance between previous separate clouds must be at least four to five grid 
intervals initially. These conditions are a combination of the definitions of merger found in 
several observational studies (Changnon 1976; Houze and Change 1977; Simpson et al. 1980). 
Based on the GCE model results, the most unfavorable environmental conditions for cloud 
merging are 1) less unstable stratification of the atmosphere and 2) weaker large-scale forcing. 

One advantage of the model simulations is that the model can be rerun in order to 
investigate the sensitivity of its results to various physical processes. For example, a pair of 
runs using identical initial conditions were performed. The only difference is that the drag 
force of rain water in the vertical equation is set to zero in the sensitivity test. The absence of 
the drag force can lead to a delay in either the onset or the weakening of the downdraft below 
the cloud. The new convective cell in the merged situation did not occur in the run with weaker 
downdrafts. This sensitivity test demonstrated the importance of downdrafts on merger. 

Later, a total of nine three-dimensional experiments were made using the same GATE 
data set (Tao and Simpson 1989a). Ten merged systems involving precipitating clouds were 
identified. Eight of these ten mergers involved two previously separated clouds (cells E and 
F); seven of these lie along a line roughly parallel to the initial environmental wind shear vector 
(called parallel cells, see Fig. 1). Only one merger lies along a line roughly perpendiculai to 
the wind shear vector prior to the merging (called perpendicular cells, see Fig. 2). The 
dominance of parallel cells is consistent with observations in FACE and GATE (Simpson et al. 
1980; Turpeinen 1982). The remaining two systems involve three clouds and are a 
combination merger of parallel and perpendicular cells. It was also found that a cloud bridge, 
which consists of a few low-level cumuli which develop and connect the clouds before the 
merger is detected on radar, occurs in most of the simulated merger cases. (This phenomenon 
was also well-simulated in the 2-D model.) New cell (Cell G in the parallel merger case and 
Cell K in the perpendicular meager case) at the cloud bridge area developed vigorously. Both 
backward and forward air parcel trajectory analyses (Fig. 3) were performed. Forward air 
parcel trajectories are computed using grid points located in the merging area. Then, a 
backward trajectory calculation was performed to locate the origins of the high-rising parcels. 
These trajectory analyses show that the high-rising air parcels at the bridge area originated close 
to or within the regions occupied by previous separated cells (Cells E and F). These air parcels 
were strongly affected by either one or two interacting cold outflows. Both 2D and 3D GCE 
model studies clearly suggest that the primary initiating mechanism for the occurrence of a 



precipitating cloud merger is the cloud downdrafts and their associated cold outflows as 
proposed by Simpson (1980). A significant difference between the simulated parallel and 
perpendicular cells is that the latter cells are usually situated closer to each other (5-6 km) prior 
to merging, compared to the former (10 km or more). An explanation for this difference is that 
the direction of individual cell movement as well as the direction of cold outflow are 
predominantly directed down shear. 

3.1.2 Comparison with Other Cloud Resolving Model Results 

The causes of merging have been investigated by Hill (1974), Wilkins etal. (1976), Orville et 
al. (1980), Turpeinen (1982), Bennetts et al. (1982) and Kogan and Shapiro (1996) using 
cloud resolving models. Orville etal. (1980) investigated the effects of varying the spacing, 
timing and intensity of two initial impulses in the context of a two-dimensional cloud model 
including warm rain and hail processes. Merging was found to result if two clouds were 
relatively close to each other (less than 7 km) and if the clouds were of different strength or 
initiated at different times (at intervals of 6 minutes). The mechanism of merging was 
attributed to the existence of a pressure gradient directed from the weaker and younger cell 
toward the older and stronger one. By using a three-dimensional cloud model, Turpeinen 
(1982) also found that the mechanism of merging was dependent on the perturbation pressure 
distribution. Note that these two modeling studies used the joining of the 100% relative 
humidity isopleth of water vapor as a criterion for merger. The formation of a cloud bridge 
observed by Simpson (1980) has been simulated by both studies. But, vigorous development 
of the new convective cell at the cloud bridge area did not occur in Orville et al. (1980) and 
Turpeinen (1982). Turpeinen (1982) suggested that this discrepancy might be attributed to the 
absence of mesoscale convergence in the model simulations. 

Kogan and Shapiro (1996) performed three-dimensional numerical simulations of 
mergers using explicit microphysics in a shear-free environment. Their criterion for cloud 
merger was based on the visual form of cloud updraft merger on a horizontal cross section. An 
arbitrary contour interval specified in the graphics routine (2 m/s) for coalescence of vertical 
velocity was used. This criterion was examined every 300 s. Kogan and Shapiro (1996) 
found that updraft merger occurred in four of the six simulations. They also found that after 
updraft merger, the maximum vertical velocity and domain averaged kinetic energy were 
increased over the single bubble simulation. They hypothesize that the mergers were a 
consequence of mutual advection, that is each of the clouds advected its neighbor in its radial 
inflow. Bennetts et al. ( 1 982) also attributed merging in their numerical simulations to "mutual 



attraction". Kogan and Shapiro (1996) also found that the most favorable conditions for 
merger occur when the cells are closer than 4.5-6 radii apart (about 3-5 km between the centers 
of the temperature perturbations). No vigorous development occurred after the two updrafts 
merged, however. No precipitating downdraft was present in their simulations. 

There is one major difference between the GCE model simulated mergers and those 
from others (Orville et al. 1980; Turpeinen 1982; Bennetts et al. 1982; Kogan and Shapiro 
1996). The simulated mergers from other modeling studies are the consolidation of two initial 
independent single bubbles (the first-order merger). Their simulated mergers do not have 
vigorous development in contrast to the GCE model simulations. The basic design of the GCE 
modeling study is to generate several convective clouds randomly inside the model domain 
and, then, to observe and analyze the interactions between the simulated clouds. Neither 
locations nor intensities of simulated clouds are predetermined. The mergers identified in Tao 
and Simpson (1984, 1989a) only involve precipitating clouds (by definition). Their merged 
cases lasted longer and produced quite significant surface precipitation as observed by Simpson 
et al. (1980, 1993). Tao and Simpson (1989a) found that some of the previously distinct 
clouds associated with merger cases resulted from the consolidation of smaller-sized clouds. 
[This may also explain why the mergers discussed in Tao and Simpson (1989a) are very 
similar to the second-order merged systems observed by Simpson et al. (1980).] These 
smaller sized cells were predominantly oriented along the direction of the wind shear vector 
when they merged together. This result is inconsistent with the simulation performed by 
Turpeinen (1982). Situations for this type of merger only involve shallow clouds with little or 
no surface precipitation. Thus, the mechanism responsible for their merging can not be cloud 
downdrafts and their associated cold outflows. The pressure distribution, as suggested by 
Orville et al. (1980) and Turpeinen (1982), mutual advection, as suggested by Kogan and 
Shapiro (1996) and the differential motions between convective elements (LeMone 1989) are 
probably the major mechanisms for this type of merger. All first-order simulated mergers may 
require is for two initially separated convective cells to be very close [from 7 km in Orville et 
al. (1982) to about 4 km in Kogan and Shapiro (1996) and Tao and Simpson (1989a)]. 

The definition of cloud merger is not unique in observational studies (Westcott, 1984). 
Observational studies are based on radar derived information. The observational studies 
usually define merger in terms of coalescence of precipitation areas or radar reflectivity (at 1 
mm h"l, minimum detectable reflectivity signal). Additional criteria related to the distance 
between initially distinct convective elements and the duration of precipitation are also 
sometimes applied. Numerical simulations have used modeled dynamical and 


thermodynamical parameters (i.e., overlap of buoyancy, updraft, humidity, or circulation 
fields) to define mergers. Westcott (1984) pointed out that in order to perform better merger 
studies, it is necessary to clearly relate convective system's dynamical, thermodynamic and 
microphysical structures and their radar image. 

3.2 Convective - Stratiform Interaction 

One of the major findings from GATE was the important contribution to rainfall from 
mesoscale convective systems (MCSs 2 ). For example, Houze (1977) estimated that four MCSs 
accounted for 50% of the rainfall at one of the GATE ships during Phase III. It was also 
estimated that the widespread stratiform rain accounted for about 32%-49% of the total rainfall 
from the GATE MCSs (Houze 1977; Zipser et al. 1981; Gamache and Houze 1983). In 
addition, observations indicated that little stratiform rain fell during the early stages of tropical 
MCSs. As the stratiform cloud developed and expanded, the total amount of rain falling from 
it became equal to that generated in the convective region. The fraction of stratiform rainfall 
from midlatitude squall lines has been estimated at 29%-43% (Rutledge and Houze 1987; 
Johnson and Hamilton 1988). The existence of unsaturated warm mesoscale descent beneath 
the stratiform region was identified by Zipser (1969) and conceptualized in Houze (1977) and 
Zipser (1977). The associated mesoscale ascent at the middle and upper layers of the stratiform 
region was diagnosed from indirect observations by Gamache and Houze (1983). One type of 
MCS is a squall line. The conceptual model of tropical and midlatitude squall lines are shown 
in Fig. 4. 

The vertical distribution of heating in the stratiform region of MCSs is also 
considerably different from the vertical profile of heating in the convective region (Houze 1982; 
Johnson 1984). The convective profiles always show heating throughout the depth of the 
troposphere which is maximized in the lowest 2-5 km. The shapes of the heating profiles are 
quite similar with only slight variations in their magnitude for different MCSs from different 
geographic locations. The same can generally be said about the stratiform region. Heating is 
maximized in the upper troposphere, however, between 5 and 9 km while cooling prevails at 
about 4 km. In addition, many recent studies (Adler and Negri 1988; Tao et al. 1993b) 
indicated that a separation of convective and stratiform clouds is necessary for a successful 
surface rain and latent heating profile retrieval from remote sensors. 


Houze (1997) defined a mesoscale convective system (MCS) as "a cloud system that occurs in 
connection with an ensemble of thunderstorms and produces a contiguous precipitation area ~ 100 km 
or more in horizontal scale in at least one direction". 
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These findings lead to an important question: what are the origins and growth 
mechanisms of particles in stratiform precipitation? Chen and Zipser (1982) suggested that 
both depositional growth associated with upward motion in the anvil and the horizontal flux of 
hydrometeors from the convective region are important in the maintenance of anvil 
precipitation. In a kinematic model study of a GATE squall line, Gamache and Houze (1983) 
showed quantitatively that 25-40 percent of the stratiform condensate was created by mesoscale 
ascent at mid-to-upper levels in the stratiform region itself. Gallus and Johnson (1991) found 
that the contribution to surface rainfall from condensation in the mesoscale updraft was 
comparable in magnitude to the transport of condensate rearward from the convective line 
during a rapidly weakening stage of a mid- latitude squall line. Using a kinematic (steady-state) 
cloud model, Rutledge (1986) suggested that the condensate produced by mesoscale ascent is 
largely responsible for the large horizontal extent of light stratiform precipitation to the rear of 
the same GATE squall line analyzed by Gamache and Houze (1983). Using higher resolution, 
Doppler-derived air motions associated with a midlatitude squall line as input in their two- 
dimensional kinematic model, Rutledge and Houze (1987) found that deposition in the 
mesoscale updraft accounted for 80 percent of the stratiform precipitation. They also 
conducted a series of sensitivity tests and found that almost no rain reached the surface in the 
stratiform region without the influx of hydrometeors from the convective cells, while only 
about one-fourth as much stratiform rain reached the surface in the absence of mesoscale 
ascent. 

3.2. 1 The GCE Model Simulations Results 

Observational studies have had to use a steady state assumption to estimate the transfer of 
hydrometeors from the convective region to its associated stratifrom region as well as a 
relatively simple 1-D cloud model to estimate the microphysical processes within the convective 
and stratiform regions. The time-dependent cloud resolving models (Tao et al. 1993a; Chin 
1994; and Caniaux et al. 1994; Tao 1995; and others) have been used to explicitly quantify the 
origins and growth mechanisms of particles in stratiform precipitation by calculating the water 
budgets (microphysical processes and transfer processes of hydrometeors between convective 
and stratiform regions). 

Several organized convective systems (EMEX, TOGA COARE, TAMEX and 
PRESTORM), which occurred in different large-scale environments, have been simulated by 
the GCE model and the associated water budgets were analyzed (Tao et al. 1993a; Tao 1995). 



Table 5 compares several characteristics of the large-scale flow (i.e., stability, Richardson 
number, and precipitable water) in which these convective systems were embedded. The 
propagation speed of these systems and the references for the GCE model simulations are also 
listed. The convective available potential energy (CAPE) associated with the tropical 
convective systems is moderate (from 1400 to 1660 m2 s'^) and smaller than that of the 
midlatitude system (PRESTORM). The vertical integrated water vapor contents are much 
higher for the TOGA CO ARE and EMEX cases compared to the PRESTORM case. 

The water budgets in the convective, stratiform, and nonraining regions associated with 
the TOGA COARE, EMEX, TAMEX and PRESTORM convective systems are shown in Fig. 
5. The water budgets are separated into three different layers: lower (surface to 10 °K level), 
middle (from 10 K to -10 K) and upper (-10 K to 100 mb). The horizontal transfer of 
hydrometeors from the convective to the stratiform region occurs mainly in the middle 
troposphere for the EMEX and TOGA COARE convective systems. By contrast, two thirds of 
the horizontal transfer of hydrometeors is accomplished in the upper troposphere for the 
PRESTORM case. This is caused by the strong convective updrafts associated with the 
PRESTORM case. Also a more vigorous transfer of hydrometeors in the lower troposphere 
from the stratiform region back into the convective region occurs for the PRESTORM case. 
This is a consequence of the strong rear inflow simulated for this midlatitude case. For the 
TAMEX case, the horizontal transfer of hydrometeors can occur in both the middle and upper 
troposphere. A downward transfer of hydrometeors from the middle to the lower troposphere 
is a dominant process in the stratiform regions for all four cases. The interaction between the 
stratiform and non-surface-raining region is less significant than that between the convective 
and stratiform region. 

The contribution to stratiform rain by the convective region has to be quantified by 
estimating a ratio (R), R=C T /(C T +C m ), where C x is the horizontal transfer of hydrometeors 
from the convective region into the stratiform region above the 10 °C level, and C m is the sum 
of the net condensation in the stratiform region and in the non-raining region above the 10 °C 
level. A small ratio indicates that the horizontal transfer of hydrometeors from the convective 
region is a small source of condensate for the stratiform anvil, whereas a ratio near unity 
indicates that nearly all of the condensate in the stratiform region was transported from the 
convective region. All four GCE modeled cases showed large ratios, from 0.33 to 0.8^., 
implying the role of the convective region in the generation of stratiform rainfall can not be 
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neglected 3 (Table 6). The relative importance of the horizontal transfer processes to the 
stratiform water budget is similar between the initial and the mature stages of the TAMEX, 
TOGA COARE and EMEX systems, and this is likely due to the fact stratiform clouds 
developed rapidly. In contrast, during the initial stage of the PRESTORM simulation, nearly 
all of the condensate in the stratiform region was a result of the horizontal transport from the 
convective region. As the PRESTORM system matured, the contribution made by the 
horizontal transport of hydrometeors from the convective region (i.e., the ratio R) decreased, 
such that the sources of condensate in the stratiform water budget were similar for all of the 
mature storms. It is hypothesized that during the initial stage of the PRESTORM simulation, 
much of the condensate transported from the convective region is used to moisten and modify 
the dry environment at middle and upper levels. Condensation and deposition become 
increasingly more important with time in the stratiform water budget once the larger-scale 
environment reaches saturation. This evolution in the stratiform water budget is less obvious 
in the TAMEX, TOGA COARE and EMEX cases because the environment is much more 
moist. (Note that the TAMEX, EMEX and TOGA COARE cases have more stratiform rainfall 
than its PRESTORM counterpart.) The GCE model results also indicated that the similarity in 
R at the mature stages of all systems is likely to be result of large stratiform regions being 
present. 

3.2.2 Comparison with Observational and Other Cloud Resolving Model Results 

Table 7 lists the ratio (R) from observational studies using composite wind and thermodynamic 
fields for five different GATE MCSs (Leary and Houze 1980; Gamache and Houze 1983), a 
midlatitude squall line (Gallus and Johnson 1991), and a tropical-continental squall line (Chong 
and Hauser 1989). For six out of the seven observed cases the ratio is very close to or above 
0.50. This implies that the convective region plays a very important role in the generation of 
stratiform rain. Very good agreement is evident between the ratio at the mature stage of the 
modeled PRE-STORM squall system and that estimated by Gallus and Johnson (1991). The 
modeled EMEX and TOGA COARE cases indicate a relatively small contribution (0.37) to 
stratiform formation from the convective region compared to those determined from the 
kinematic studies. 

3 The convective region can also transport the water vapor originally from low troposphere 

into the stratiform region. Sui et al. (1994) indicated that this water vapor transport is the source 
for stratiform formation (deposition and condensation). Dynamic triggering of the stratiform 
formation can be gravity waves excited by strong and deep convective cells associated with the 
convective region. 



Table 7 also shows the ratios determined from other CRM results (Chin 1994; Chin et 
al. 1995; Caniaux et al. 1994). Very good agreement is evident between the ratio at the mature 
stage of our modeled EMEX and TOGA COARE squall systems and a tropical squall case 
simulated by Chin et al. (1995). The comparison between our simulated PRE-STORM and 
other CRM simulated midlatitude cases (Chin 1994; and Caniaux et al. 1994), however, is 
quite different. Caniaux et al. (1994) suggested that the smaller contribution from the 
convective region to stratiform formation compared to observational studies (Chong and 
Hauser 1989) was due to the inability of convective updrafts to transport condensate to high 
levels in their two-dimensional simulation, the slower propagation speed and the existence of a 
transition zone. The very smaller stratiform portion (10%) in the midlatitude case simulated by 
Chin (1994) is the reason for the higher R. The GCE modeled PRESTORM case has a higher 
R (0.8) and smaller stratiform portion (14%) at the initial stage. 

A direct comparison between these studies and the GCE model studies should be done 
with caution, because a different spatial resolution and a different definition for the convective- 
stratiform regions was used. For example, Caniaux et al. (1994) had a fixed number of model 
grid points (50) designated as the convective region. Remaining grid points with surface 
precipitation comprised the stratiform region. In the GCE model, the convective and stratiform 
regions are identified using information from surface rainrates first (i.e., Churchill and Houze 
1984). Additional criteria are applied which have been included to identify regions where 
convection may be quite active aloft though there is little or no precipitation yet at the surface, 
such as areas associated with tilted updrafts and new cells initiated ahead of organized squall 
lines (Tao et al. 1993, Lang et al. 2000). The GCE method was adopted by Chin (1984) and 
Chin et al. (1995). The comparison between the GCE simulated PRE-STORM, TOGA 
COARE, TAMEX and EMEX cases, however, is consistent because the same type of data set 
and the same criteria for partitioning the convective and stratiform regions were being used. 

3.2.3 The Convective and Stratiform Processes in Large-Scale Models 

Molinari and Dudek (1992) and Frank (1993) suggested that the best approach to cumulus 
parameterization in large-scale scale models (30-120 km horizontal resolution, 150-300 second 
time steps) appears to be "to use a scheme that operates simultaneously with and interacts 
explicitly with the explicit scheme (grid scale microphysical processes)". They termed such 
schemes " hybrid schemes". The hybrid approach (by separating out the forcing mechanism for 
the mesoscale component) resolves the "mesoscale" circulations and microphysical processes 



15 


that directly influence the development of the "stratiform clouds". Cumulus parameterization 
makes use of steady-state cloud models that interact with grid-scale variables and provide net 
heating, drying, and condensate associated with "convective cells". The interaction between 
parameterized and explicitly resolved cloud processes is through the detrainment of water vapor 
and condensate generated from the steady-state cloud model into the "resolved" stratiform 
clouds. 


Recently, GCMs and climate models (i.e., CSU GCM and GISS GCM) allow both a 
cumulus parameterization scheme and an explicit moisture scheme to be activated 
simultaneously in the model simulations. The cumulus parameterization scheme is generally 
used to represent convective precipitation (10 km spatial scale) and the explicit moisture scheme 
to represent grid-resolvable precipitation such as stratiform/cirrus clouds (100-200 km spatial 
scale). The CSU GCM has implemented an explicit microphyiscal scheme with five 
prognostic variables for the mass of water vapor, cloud water, cloud ice, rain and snow 
(Fowler el al. 1996). The GISS global climate model has added an efficient prognostic cloud 
water (one species only). Stratiform clouds can be coupled with parameterized convection 
through detrainment of cloud water and/or cloud ice from the "tops" of cumulus towers or at 
any level above 550 mb (Del Genio et al. 1996). 

The explicit interaction between cumulus parameterized and grid-scale resolved 
microphysics is only one-way in the current large-scale models. Note that some water 
condensate generated by the stratiform region can be transported into the convective region. In 
addition, how much (all or partial) and where the (cloud tops or above melting layer) 
parameterized water condensate should detrain into the explicitly resolved microphysical 
scheme needs to be addressed. The CRM results can and should be used for improving the 
cumulus parameterization schemes as well as for understanding the interaction between the 
cumulus parameterization schemes and the explicit moisture schemes. In the future, CRMs can 
be used to study the time evolution of each of the water budget terms associated with MCSs in 
different geographical regions as well as to determine whether any important variations in the 
evolution of the water budget can be explained in terms of differences in the wind and 
thermodynamic characteristics of the large-scale environments. 

3.3 Mechanisms of Cloud-Radiation Interaction 

The interaction between clouds and radiation is two-way. On the one hand, clouds can reflect 
incoming solar and outgoing long-wave radiation. On the other hand, radiation can enhance or 
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reduce the cloud activity. Gray and Jacobsen (1977) suggested that differential cooling between 
cloudy and clear regions can enhance cloud activity in the cloudy region. Long-wave radiation 
cools the stratiform cloud top but warms the stratiform cloud base (Cox and Griffith 1979). As 
a result, long-wave radiation can destabilize the stratiform cloud layer. Webster and Stephens 
(1980) also suggested that this destabilization was quite an important process in the light 
precipitation region during WMONEX. Stephens (1983) further suggested that the effects of 
radiation on the growth and sublimation rates of ice particles are significant. Particle growth 
(sublimation) is enhanced (suppressed) in a radiatively cooled (heated) environment. Radiative 
cooling could also destabilize the large-scale environment (Dudhia 1989). The cloud-radiation 
interaction can also have major impact on the diurnal variation of precipitation processes over 
the tropics. For example, the thermodynamic response of clouds to radiative heating [cloud 
development is reduced by solar heating and enhanced by IR cooling - Kraus (1963), Randall 
et al. (1991)] and the large-scale dynamic response to the radiational differences between 
cloudy and clear regions (Gray and Jacobson 1977) have been suggested as the mechanisms 
responsible for the diurnal variation of precipitation over tropical oceans. 

3.3.1 The GCE Model Simulated Results 

A two-dimensional version of the GCE Model has been used to perform a series of sensitivity 
tests to identify which is the dominant cloud-radiative forcing mechanism with respect to the 
organization, structure and precipitation processes for both a tropical (EMEX) and a midlatitude 
(PRESTORM) mesoscale convective system (Tao et al. 1996). Figure 6 shows the schematic 
diagram demonstrating the impact of cloud-radiation mechanisms on surface precipitation for 
EMEX and PRESTORM cases. The GCE model results indicated that the dominant process 
for enhancing the surface precipitation in both the PRE-STORM and EMEX squall cases was 
the large-scale radiative cooling. However, the overall effect is really to increase the relative 
humidity and not the CAPE. Because of the high moisture in the tropics, the increase in 
relative humidity by radiative cooling can have more of an impact on precipitation in the tropical 
case than in the midlatitude case. The large-scale cooling led to a 36% increase in rainfall for 
the tropical case. The midlatitude squall line with a higher CAPE and lower humidity 
environment was only slightly affected (7%) by any of the longwave mechanisms. The GCE 
model results also indicated that the squall systems' overall (convective and stratiform) 
precipitation is reduced by turning off the cloud-top cooling and cloud-base warming. 
Therefore, the cloud-top cooling - cloud-base warming mechanism was not the responsible 
cloud-radiative mechanism for enhancing the surface precipitation. However, the circulation as 
well as the microphysical processes were indeed (slightly) enhanced in the stratiform region by 
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the cloud-top cooling and cloud-base warming mechanism for the midlatitude squall case. For 
both cases, the model results show that the mechanism associated with differential cooling 
between the clear and cloudy regions may or may not enhance precipitation processes (+5% to 
-7%, respectively for the EMEX and PRESTORM cases). However, this mechanism is 
definitely less important than the large-scale longwave radiative cooling. 

Solar heating was run from 9 AM to 1 PM LST in both environments and was found to 
decrease the precipitation by 7% in each case, compared to the runs with longwave radiation 
only. This result suggests that solar heating may play a significant role in the daytime 
minimum/nighttime maximum precipitation cycle found over most oceans, as noted in the 
observational study of Kraus (1963). Sui et al. (1998) used the GCE model and performed a 
15 day integration to simulate TOGA COARE convective systems. Their simulated diurnal 
variation of surface rainfall is in reasonable agreement with that determined from radar 
observations. They also found that modulation of convection by the diurnal change in available 
water as a function of temperature was responsible for a maximum in rainfall after midnight. 
This simply implies that the increase (decrease) in surface precipitation associated with IR 
cooling (solar heating) was mainly due to an increase (decrease) in relative humidity. The GCE 
model results also showed that the diurnal variation of sea surface temperature only plays a 
secondary role in diurnal variation in precipitation processes. 

3.3.2 Comparison with Other Cloud Resolving Model Results 

Table 8 lists the previous modeling studies that have investigated the impact of cloud-radiation 
interactive processes on various cloud systems. The increments in surface precipitation in 
Table 8 are normalized against the run without radiative processes. The conclusions associated 
with cloud-radiation mechanisms for our GCE modeled tropical (EMEX) and mid-latitude 
(PRESTORM) squall cases are in good agreement with many of these previous modeling 
studies. For example, Xu and Randall (1995), Miller and Frank (1993) and Fu et al. (1995) 
indicated that the differential cooling between cloudy and clear regions plays only a secondary 
role in enhancing precipitation processes. Xu and Randall (1995) and Fu et al. (1995) 
suggested that the cloud-top cooling and cloud-base warming destabilization mechanism could 
belmportant for prolonging the lifespan of high anvil clouds (around 10 km). Xu and Randall 
(1995) showed that this direct cloud destlbilization does not have any impact on surface 
precipitation. The modeling studies (Fu et al. 1995; Miller and Frank 1993) also indicated that 
more surface precipitation can be generated in runs with constant clear-air radiative cooling than 
without. In addition, previous modeling results (Chin 1994; Chin et al. 1995; Miller and Frank 
1993) indicated that solar radiative processes can reduce precipitation processes. However, the 
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amount of increase or decrease in surface precipitation varies quite significantly among these 
different modeling studies, but only in regard to the tropical convective systems and not the 
midlatitude systems. One possible explanation is that large-scale forcing (lifting) was needed 
in some of these different tropical convective system studies. The imposed lifting varied from 
2 cm/s to 14 cm/s in magnitude and was applied continuously or discontinuously in time 
among the different studies (see Table 8). Using an earlier version of the GCE model (Tao and 
Simpson 1989b), which included a superimposed large-scale vertical velocity as the main 
forcing, sensitivity tests using two different large-scale vertical velocities were performed. The 
results show that the radiative effects on the clouds are quite sensitive to the imposed 
background ascent (or lifting). The larger the imposed vertical velocity (9-12 cm/s) there is, 
the lesser the impact of longwave cooling on surface precipitation processes (over 24 h of 
simulation time). Frank and Miller (1993) also obtained a similar conclusion using a regional 
scale model. Also note that the larger the imposed vertical velocity, the larger the cloud 
coverage that was generated. 

The physical processes responsible for diurnal precipitation were found to be quite 
different between the GCE model and other CRM studies. For example, Xu and Randall 
(1995) found that nocturnal convection has basically a direct result of cloud-radiation 
interactions, in which solar absorption by clouds stabilized the atmosphere. Their simulated 
rainfall for both non-interactive and interactive radiation were quite similar, however. Liu and 
Moncrieff (1998) showed that direct interaction of radiation with organized convection was the 
major process that determined the diurnal variability of rainfall. Their results indicated that well 
(less) organized cloud systems can have strong (weak) diurnal variations of rainfall. They also 
suggested that ice processes are needed. The model set-ups between Sui et al. (1998) and Liu 
and Moncrieff (1998) are quite different, however. In Liu and Moncrieff (1998), the 
horizontal momentum was relaxed to its initial value that had a strong vertical shear in 
horizontal wind. On the other hand, the horizontal wind was nudged to time-varying observed 
values in Sui et al. (1998). Consequently, only long-lived squall lines (or fast-moving 
convective systems) were simulated in Liu and Moncrieff (1998) over the entire simulation. In 
Sui et al. (1998), however, their simulated cloud systems had many different sizes and various 
life cycles. A more rigorous cloud resolving model inter-comparison involving mechanisms 
associated with diurnal variation is needed in the future. A good quality controlled long term 
observational data set that can provide large-scale initial conditions is also required. 


3.4 Latent Heating Profiles and TRMM 
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The GCE model has been used to develop a Convective-Stratiform Heating (CSH) algorithm. 
The CSH algorithm uses surface precipitation rates, amount of stratiform rain, and information 
on the type and location of observed cloud systems as input. The CSH algorithm also utilizes a 
lookup table that consists of convective and stratiform diabatic heating profiles for various 
types of cloud systems in different geographic locations. These profiles are obtained from 
GCE model simulations by temporally and spatially averaging the heating distributions in the 
convective and stratiform regions of the systems, which are then normalized by their total 
surface rainfall (see Fig. 7). The heating profiles (normalized with surface rainrate) shown in 
Fig. 7 all have a characteristic shape for the convective and stratiform heating (e.g., Houze, 
1982; 1997). These include maximum convective heating in the lower to middle troposphere, 
maximum stratiform (anvil) heating in the upper troposphere, and regions of stratiform cooling 
prevailing in the iower troposphere. Also, larger heating aloft in the stratiform region is 
associated with larger cooling in the iower troposphere. However, some notable differences 
do exist. For example, the level separating the heating and cooling in the stratiform region 
(indicating the freezing or melting level) is different for the convective systems simulated by the 
GCE model and determined by diagnostic budget. The differences in the height of the 
strafiform region cooling probably reflect differences in melting layer height or the type of 
convective systems (system has erect updrafts has higher height). The cooling is quite strong 
near the surface for the African convective system due to a dry boundary layer (Caniaux et al. 
1994). The latent heating profiles modeled by the GCE and determined kinematically 4 are quite 
different for the GATE convective system. Nevertheless, there is, perhaps, more similarity 
than difference in these profiles shown in Fig. 7. This may imply that the look-up table may 
not need a significant number of heating profiles. 

Tao et al. (2000a) evaluated the CSH algorithm's performance by retrieving the latent 
heating profiles associated with three TOGA COARE convective episodes (December 10-17 
1992; December 19-27 1992; and February 9-13 1993). The inputs for the CSH algorithm 
were SSM/I (similar to TMI) and Radar (similar to TRMM PR) derived rainfall and stratiform 
amount. Diagnostically determined latent heating profiles calculated using 6 hourly soundings 
were used for validation. The temporal variability of retrieved latent heating profiles using 
radar estimated rainfall and stratiform amount was in good agreement with that diagnostically 
determined for all three periods. However, less rainfall and a smaller stratiform percentage 
estimated by radar resulted in weaker (underestimated) latent heating profiles and lower 


4 The convective and stratiform heating profiles were derived using composite "kinematic 

and thermodynamic" fields from radar, upper air soundings and aircraft measured winds. 
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maximum latent heating levels compared to those determined diagnostically. Rainfall 
information from SSM/I can not retrieve individual convective events due to poor temporal 
sampling. Nevertheless, this study suggested that a good rainfall retrieval from SSM/I for a 
convective event can lead to a good latent heating retrieval. 

The four-dimensional latent heating structure over the global tropics for February 1998 
was obtained using TRMM rain products in Tao et al. (2000b). Figure 8 shows monthly 
(February 1998) mean latent heating at three different altitudes (2, 5 and 8 km) over the global 
tropics from the CSH algorithm. The horizontal distributions or patterns of latent heat release 
identify the areas of major convective activity (i.e., a well defined ITCZ in the Pacific, a 
distinct SPCZ) in the global tropics. A well defined ITCZ in the east and central Pacific and in 
the Atlantic Ocean, a distinct S. Pacific Convergence Zone (SPCZ), and broad areas of 
precipitation events spread over the continental regions, are present. Also, stronger latent heat 
release (10 K/day or greater) in the middle and upper troposphere is always associated with 
heavier surface precipitation. Heating in the upper troposphere over the Pacific and Indian 
Oceans is much stronger than the heating over Africa, S. America and the Atlantic Ocean. The 
difference in retrieved convective and stratiform properties between the various geographic 
locations is the major reason for the difference in the heights of the maximum latent heating 
level. Higher stratiform amounts always contribute to higher maximum latent heating levels. 
Whether the higher stratiform proportions and more frequent vigorous convective events in the 
Pacific are related to the warmer SSTs needs to be studied using multi-season and multi-year 
retrieved latent heating profiles. Note that differential heating between land and ocean in the 
upper troposphere could generate strong horizontal gradients in the thermodynamic fields and 
interact with the global circulation. 

One interesting result from Fig. 8 is the relatively strong cooling (-1 to -2 K/day) at 2 
km over the (East, Central and South) Pacific and Indian Oceans but not over the continental 
regions (i.e., Africa and S. America). This result is due to the fact that the TMI observations 
had less stratiform precipitation over the continental regions which is not conducive to 
retrieving stronger low level cooling over the continental regions relative to the tropical oceans. 
However, it is still not an expected result because the moisture content is higher over oceans. 
Cooling by evaporation of raindrops in the lower troposphere should be stronger over dry 
areas. Several previous observational studies were performed to analyze the heating budget 
obtained from sounding networks over the Pacific warm pool region and the Amazon region. 
For example, Lin and Johnson (1996) found weak cooling at low levels, probably induced by 
mesoscale downdrafts or evaporation by shallow cumuli, in the mean heating profile over the 
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TOGA CO ARE region for the month of February 1993. Greco et al. (1994) calculated latent 
heating profiles from the ABLE network. Their results indicated that the distribution of heating 
is quite similar to the studies of those of West African squall lines (Chong and Hauser 1990). 
Peak heating occurs between 500 and 550 hPa (about 5-6 km). Their results did not exhibit 
low level diabatic cooling for the ABLE case. They suggested that the lowermost 2-3 km over 
the Amazon rain forest canopy is characterized by a strong diurnal cycle of evapotranspiration 
and upward convective fluxes of moisture producing very large mixing ratios (Fitzjarrald et al. 
1990). Model results (Scala et al. 1990) also suggested that dry tropospheric air is not present 
for the production and maintenance of evaporatively cooled downdrafts. The high moisture 
content during the wet season in the lower troposphere of the Amazon Basin may prevent or 
severely limit cooling below cloud base. Thus, more low level cooling over the Pacific than 
over S. America as estimated by the CSH heating algorithm is, perhaps, reasonable. 

The CSH algorithm estimated heating profiles show one maximum heating level, and 
the level varies between convective activity from various geographic locations. These features 
are in good agreement with the heating profiles obtained from the results of diagnostic studies 
over a broad range of geographic locations (Yanai et al. 1973; Johnson 1984, 1992; Thompson 
etal. 1979; Houze 1989; Frank and McBride 1989; Greco et al. 1994; Frank et al. 1996; Lin 
and Johnson 1996 and many others). The magnitude of their estimated latent heating release is 
also in good agreement with those determined from diagnostic budget studies. 

Two other latent heating retrieval algorithms, the Goddard Profiling (GPROF) heating, 
and the Hydrometeor heating (HH), were also used to estimate the latent heating for February 
1998 and their results were compared to the those estimated by the CSH algorithm. The 
horizontal distribution or patterns of latent heat release from the three different heating retrieval 
methods are quite similar. They all can identify the areas of major convective activity [i.e., a 
well defined Intertropical Convergence Zone (ITCZ) in the Pacific, a distinct SPCZ] in the 
global tropics. The magnitude of their estimated latent heating release is also in good 
agreement with each other and with those determined from diagnostic budget studies. 
However, the major difference among these three heating retrieval algorithms is the altitude of 
the maximum heating level. The CSH algorithm estimated heating profiles only show one 
maximum heating level, and the level varies between convective activity from various 
geographic locations. These features are in good agreement with diagnostic budget studies. 
By contrast, a broader maximum of heating, often with two embedded peaks, is generally 
derived from applications of the GPROF heating and HH algorithms, and the response of the 
heating profiles to convective activity is less pronounced. Also, GPROF and HH generally 
yield heating profiles with a maximum at somewhat lower altitudes than CSH. 
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Heating profiles for the TRMM Field Campaign sites (i.e., SCSMEX, May - June 
1998; LB A - TRMM/Brazil, January - February 1999; and KWAJEX, July - September 1999) 
as well as other major field campaigns such as DOE/ ARM will be produced using the three 
different heating algorithms, and these will be compared to profiles determined from the field 
campaign sounding networks. This future comparison can provide an assessment of the 
absolute and relative errors of the heating retrieval algorithms. In addition, global analyses will 
be used to identify/compare the large-scale circulation patterns for the retrieved periods and for 
periods during previous field campaigns (i.e., TOGA COARE and GATE). It is reasonable to 
assume that the latent heating structures for Westerly Wind Bursts (WWBs) and Super Cloud 
Clusters (SCCs) occurring in similar large-scale circulations and with similar SSTs may not be 
very different. 

3. 5 Response of Tropical Deep Cloud Systems to Large-Scale Processes 

Generally speaking, the role of clouds in the atmospheric general circulation and global climate 
is two-fold. On the one hand, clouds owe their origin to large-scale dynamical forcing, 
radiative cooling in the atmosphere and turbulent transfer processes between the ground and the 
atmosphere (e.g., the transfer of heat and moisture from the ocean to the atmosphere). On the 
other hand, the latent heat from precipitating clouds provides most of the energy received by 
the atmosphere. Clouds also serve as important mechanisms for the vertical re-distribution of 
momentum, trace gases (including Greenhouse gas, CO 2 ) and sensible and latent heat on the 
large-scale. They also influence the coupling between the atmosphere and the earth's surface 
as well as the radiative and dynamical-hydrological balance. 

The use of cloud resolving models (CRMs) in the study of tropical convection and its 
relation to the large-scale environment can be generally categorized into two groups. The first 
approach is so-called " cloud ensemble modeling" . In this approach, many clouds of different 
sizes in various stages of their lifecycles can be present at any model simulation time. The 
large-scale effects which are derived from observations are imposed into the models as the 
main forcing, however. In addition, the cloud ensemble models use cyclic lateral boundary 
conditions (to avoid reflection of gravity waves) and require a large horizontal domain (to allow 
for the existence of an ensemble of clouds). The clouds simulated from this approach could be 
termed "continuous large-scale forced convection" . This approach is always applied for 
simulation associated with tropical deep convection. On the other hand, the second approach 
for cloud resolving models does require large-scale effects to initialize the clouds. This type of 
approach usually requires initial temperature and water vapor profiles which have a medium to 


large CAPE, and an open lateral boundary condition is used. The modeled clouds, then, are 
initialized with either a cool pool, warm bubble or surface processes (i.e., land/ocean fluxes). 
These modeled clouds could be termed " self-forced convection" . The key developments in the 
cloud ensemble modeling using the continuous large-scale forced convection approach over the 
past two decades were listed in Table 2. 

3.5. 1 Simulated Results from the GCE Model 

Tao et al. (1987) used 2-D and 3-D versions of the GCE model to study the statistical 
properties of cloud ensembles for a well-organized intertropical convergence zone (ITCZ) 
rainband that occurred during GATE. The statistical properties of clouds, such as mass flux by 
cloud drafts and vertical velocity as well as condensation and evaporation were examined. 
Figure 9 shows heating rates by condensation (c) and evaporation (e) in the 2-D model and the 
3-D model. The heating rate estimated from large-scale observations, Qj 5 - Qr, and the total 
cloud heating rate in the 3-D model are also included. The rate of condensation in the 3-D case 
is slightly larger than its 2-D counterpart, but so is the rate of evaporation. As a consequence, 
the net heating effect of clouds in the 3-D model is nearly equal to the 2-D counterpart, and 
both agree with Ql - Qr as estimated from the large-scale heat budget. The GCE model also 
found that the 3-D modeled surface rainfall rates have smaller standard deviation than their 2-D 
counterparts. Overall, the GCE model results indicated that collective thermodynamic feedback 
effects and vertical transports of mass, sensible heat and moisture by the convective cells show 
profound similarities between the two- and three-dimensional GCE model simulations. 

Zipser and LeMone (1980) and LeMone and Zipser (1980) presented the results of 
statistical analyses of convective updrafts and downdrafts. Their analyses were based on 
aircraft data gathered from cumulonimbus cloud penetrations for six days during GATE. In 
order to facilitate a comparison between our model results and their analysis results, we 
subdivided the updrafts and downdrafts into active or inactive updrafts and downdrafts (Table 
9). For example, a grid point in the model is designated as an active updraft region if (a) the 
total liquid water content exceeds 0.01 g kg' 1 and (b) the vertical velocity is larger than 1ms' 1 
(or 2 m s' 1 , depending upon how we define "active") at that grid point and at that integration 
time. These ratios between active cloud updrafts and downdrafts indicate an excellent 
agreement among results from the 2-D and 3-D models as well as the cores as measured by 
Zipser and LeMone (1980). Both 2-D and 3-D model results also showed a similar feature in 
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Ql is the apparent heat source budget defined in Yanai et al. (1973). 



that the active updrafts account for approximately 75% of the upward mass flux due to clouds 
and yet they only cover about 12-14% of the total area. This result is consistent with the 
concept, first proposed by Riehl and Malkus (1958; see also Riehl and Simpson 1979), that hot 
towers play a critical role in the heat and moisture budgets in the tropics, even though they 
occupy a small fraction of the area. Overall, our comparison study has indicated that the 
statistical properties of the clouds obtained in the 2-D model are essentially the same as the 3-D 
counterpart given an identical large-scale environment (see Tao 1983). The explanation for this 
similarity between the 2D and 3D simulations is that the same large-scale advective forcing in 
temperature and water was superimposed into the GCE model as the main forcing. The cyclic 
lateral boundary condition used in the GCE model does not allow for additional forcing in the 
model domain. A two-dimensional simulation should, therefore, give a good approximation of 
the continuous large-scale forced convection. 

Large-scale models (i.e., general circulation and climate models) require not only the 
global surface rainfall pattern but also the associated vertical distribution within the Qi and Q2^ 
budgets. The GCE model can help to identify which processes should be parameterized by the 
large-scale model, as well as provide information on the vertical profiles of the Qi and Q 2 
budgets (Tao 1978; Soong and Tao 1980; Tao and Soong 1986; Tao et al. 1993a and many 
others). The GCE model was used to examine the Qi and Q 2 budgets of various cloud 
systems that developed in different geographic locations (GATE, EMEX, PRESTORM, TOGA 
COARE, ABLE, TAMEX and others). In all of these simulations, the heating due to the 
vertical eddy convergence/divergence term in sensible heat by convective clouds is always one 
order of magnitude smaller than that produced by condensation at most levels [Fig. 10(a)], On 
the other hand, the maximum value of the cooling rate by evaporation is more than half of the 
heating rate by condensation. This finding implies that the sum of the condensation and 
evaporation would provide a good approximation to the total cloud heating rate. The cloud 
heating effect would be considerably overestimated if heating by condensation alone is 
considered, ignoring cooling by evaporation. For the Q 2 budget, the GCE model results 
indicated that the net vertical eddy convergence/divergence of moisture by clouds is generally 
smaller than the rate of condensation or evaporation, but it is not negligible [Fig. 10(b)], The 
different roles of the vertical eddy convergence/divergence term in the Qi (temperature) and Q 2 
(water vapor) budgets is the major reason for Qi and Q 2 decoupling (the level of maximum 
values in the Qi and Q 2 profiles is not at the same level). The GCE model generated heating 
and drying effects agree well with those estimated from observations. 
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Q2 is the apparent moisture sink budget defined in Yanai et al. (1973). 
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GATE (1974 - in the East Atlantic) and TOGA CO ARE (1992-1993 - in the West 
Pacific warm pool region) are perhaps two of the best planned and well coordinated field 
campaigns for understanding tropical convective systems and their interactions with the large- 
scale environments within which they are embedded. An improved GCE model (ice 
microphysics, cloud-radiation interaction, dynamics, and surface fluxes) was used to study the 
response of convective systems to the large-scale environment using the data collected during 
GATE and TOGA COARE. The GCE model was integrated for 7 and 8 days, respectively, for 
GATE (September 1 to 8, 1974) and TOGA COARE (December 19 to 27, 1992). Both runs 
used 1024 horizontal grid points with a 1 km resolution and 41 vertical grid points with 
varying resolution (40 m near the surface to 1000 m at the top level). The time step was 7.5 s. 

The large-scale environments associated with the organized cloud systems that occurred 
in TOGA COARE and GATE were quite different. The large-scale advective forcing in 
temperature and water vapor as well as the large-scale vertical velocity are stronger for TOGA 
COARE than for GATE. The large-scale vertical velocity shows a diurnal signature in TOGA 
COARE but not in GATE. The mean CAPE is larger in GATE than in TOGA COARE. The 
SST is higher for TOGA COARE (about 29 °C vs 27.4 °C for GATE). The vertically 
integrated water vapor content (precipitable water) is much drier for GATE (2.47 g cm~2) than 
TOGA COARE (5.15 g cm‘2). The mean vertical shear from the surface to 700 mb of the 
large-scale horizontal wind is slightly larger for GATE than TOGA COARE during the GCE 
model simulation periods. However, the shear is much stronger from over the entire depth of 
the troposphere in TOGA COARE. The low-level wind shear can determine the organization 
of convective systems. 

Figures 1 1 (a) and (b) show the temporal variation of the GCE model simulated domain 
mean surface rain rate for TOGA COARE and GATE, respectively, There are more convective 
systems simulated by the GCE model for TOGA COARE than for GATE. This is due to the 
stronger large-scale forcing imposed in the TOGA COARE simulation. The model simulated 
surface precipitation showed a very complex structure for TOGA COARE compared with 
GATE. Overall, the GCE model-simulated cloud systems propagated in one direction while 
the individual cells embedded within the systems propagated in the opposite direction. In 
addition, the cloud tops propagate in the opposite direction of the associated surface 
precipitation. These two hierarchies of convection organization are in good agreement with 
other modeling studies (Wu et al. 1998) and with satellite observations (Nakazawa 1988; Sui 
and Lau 1989). In the GATE simulation, only shallow convective systems developed during 



the first day. Then, deep convective clouds and non-squall (slow moving) cloud systems 
developed and propagated westward with the mean wind. Squall line type (fast moving) cloud 
systems developed after September 4. After September 6, the systems simulated by the GCE 
model were less organized and produced less surface precipitation compared to the non-squall 
and squall systems. These GCE model simulated GATE features are in good agreement with 
other modeling studies (Grabowski et al. 1996; Xu and Randall 1996) and observations. 

The GCE model simulated domain averaged surface rainfall (mm), and stratiform amount 
(percentage) for both TOGA COARE and GATE are shown in Table 10. The ratios between 
evaporation and condensation, sublimation and deposition, and deposition and condensation 
were examined for both cases. These ratios illustrate the relative importance of warm verse ice 
processes and source and sink terms associated with water vapor over the course of the TOGA 
COARE and GATE simulations. The microphysical processes are broken down according to 
convective organization (i.e., slow-moving, fast-moving, less organized convective episodes 
from GATE, vigorous deep convection and weaker convective events during the Westerly 
Wind Burst period) in Table 3.5.2. As expected, more surface rainfall was simulated by the 
GCE model for TOGA COARE (153.9 mm) than for GATE (91.46 mm). Also, a higher 
stratiform component was simulated for TOGA COARE (45%) than for GATE (32%). The 
surface rainfall and stratiform component simulated by the GCE model for TOGA COARE are 
in reasonable agreement with the rainfall determined from soundings and the stratiform amount 
measured by radar [see Tao et al. (2000a) for a detailed comparison]. This close agreement is 
mainly caused by the fact that the GCE model was forced by large-scale tendencies in 
temperature and water vapor that were derived from the sounding network. However, the 
GCE model simulated surface rainfall is almost twice that estimated by radar. Johnson and 
Ciesielski (2000) indicated that the ship radars were located within a relatively dry region of the 
IFA. The lower rainfall estimates from the ship radars could also be caused by the specific Z-R 
relationship applied in Short et al. (1997). Based on radar observations (Houze 1997), the 
GCE model may have underestimated the stratiform rain for GATE fast-moving squall 
systems. The dominance of warm rain processes in the GATE squall and non-squall 
convective systems may explain the smaller stratiform rain amounts simulated by the GCE 
model. Very little ice processes on September 6 and 8 are an indication of shallower 
convection. In contrast, ice processes are quite important for both active and relatively inactive 
convective periods during TOGA COARE. Even though the large-scale environment is drier 
for GATE than TOGA COARE, evaporation is only 54% of the condensation in GATE 
compared to 71% in TOGA COARE. Weak convective episodes in both GATE and TOGA 
COARE had high ratios between evaporation and condensation compared to more intense 
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convective periods. The ratio of sublimation to deposition was smaller in the GATE 
simulation. 

3.5.2 Comparison with Other Cloud Resolving Model Results 

Krueger (1988) developed a two-dimensional cloud model with a third-moment turbulence 
closure for simulating an ensemble of cumulus clouds. He simulated the response of cumulus 
clouds to large-scale forcing, under large-scale conditions observed during GATE. Krueger 
(1988) found that cloud-scale vertical transport of moisture and evaporative cooling are 
significant in the Q2 and Qi budgets, respectively. The cloud scale vertical advection of heat is 
only important in the subcloud layer in the Qi budget. These results are consistent with our 
GCE model simulations. Lafore and Redelsperger (1991) applied a two-dimensional cloud 
model to simulate a fast moving Tropical squall line observed during COPT8 1 and a frontal 
system observed during European MFDP/FRONTS87. Their results also indicated the 
importance of evaporative cooling and cloud transport of moisture for these two cases. 
Furthermore, their results showed a relatively small effect by cloud transport of heat on the Qi 
budget except near the subcloud layer. The different roles of the vertical eddy 
convergence/divergence term in the Qi (temperature) and Q2 (water vapor) budgets is also the 
major reason for Qi and Q2 decoupling in both systems as indicated by Lafore and 
Redelsperger (1991). 

Grabowski et al. (1998) examined the effects of resolution and the third spatial 
dimension for cloud systems observed during Phase III of GATE (September 1 to 7, 1974). 
Xu and Randall (1996) used the two-dimensional model developed by Krueger (1988) to 
simulate cloud systems observed during Phase III of GATE (September 1 to 18, 1974). Wu et 
al. (1998) also used a two-dimensional model to examine the cloud properties associated with 
cloud systems observed during TOGA CO ARE (December 5 1992 to January 12 1993). 
Donner et al. (1999) used a three-dimensional model developed by Lipps and Helmer (1986) to 
simulate several GATE convective systems. The major difference for these modeling studies 
(and the improved GCE model simulation shown in Fig. 11) from the previous CRM 
simulations is that they performed long term integrations. All these studies' simulated Qi and 
Q2 budgets are in good agreement with observations. This is due to the fact that the observed 
large-scale advective forcing in temperature and water vapor were imposed as suggested by 
Soong and Tao (1980). Cloud organization in both studies also agreed well with observations 
due to the fact that the modeled simulated horizontal wind was relaxed to the observed time 
varying large-scale horizontal wind. The importance of vertical shear of the large-scale 



horizontal wind on the organization of tropical convective systems was recognized in 
theoretical studies (i.e., Moncrieff 1992) and numerical simulations (Tao and Soong 1986; 
Dudhia etal. 1987). 

Larger temporal variability in the two-dimensional integration than in the three- 
dimensional integration was found in Grabowski et al. (1998) and Donner et al. (1999), 
Donner et al. (1999) suggested that this is probably related to the different behavior of the 
CAPE and convective inhibition (CIN) in two and three dimensions. Grabowski et al. (1998), 
however, concluded that, as long as high-frequency temporal variability is not of primary 
importance, low-resolution two-dimensional simulations can be used as realizations of tropical 
cloud systems for addressing the climate problem and for improving and testing cloud 
parameterizations for large-scale models. This conclusion is only valid for CRMs using large- 
scale advective forcing and applied with periodic lateral boundary conditions. A similar 
conclusion was also obtained the results of the GCE model. 

However, there are several notable differences between two- and three-dimensional 
CRM simulations. For example, a weaker convective updraft and a stronger convective 
downdraft velocity were simulated for a GATE fast moving system in the GCE two- 
dimensional model compared with in the three dimensional model. Yet, the total upward and 
downward mass fluxes are almost identical between the two- and three-dimensional GCE 
model simulations. Lipps and Helmer (1986), however, found that their two-dimensional 
model had stronger upward and downward mass fluxes than their three-dimensional 
counterpart for the same GATE simulation as Tao and Soong (1986). They also found more 
evaporation of cloud water in the two dimensional simulation and consequently less cloud 
water was present. These results are very different from those of Wu and Moncrieff (1998) for 
simulations of TOGA COARE convective systems. More ice water and liquid water were 
simulated in the two-dimensional model than the three-dimensional model in Wu and Moncrieff 
(1998). The different cases simulated between Lipps and Helmer (1986) and Wu and 
Moncrieff (1998) is, perhaps, one of the major reasons for the differences. The microphysical 
schemes used in these two studies are also different. More detailed comparisons are needed in 
future. The GCSS model intercomparison project and field campaigns (ARM, TRMM LBA. 
TWAJEX and NASA CAMEX 7 ) can provide good quality observational data sets for CRM 
initialization as well as for its validation. 


7 ARM stand for Atmospheric Radiation Measurement, LBA for Large Scale 
Biosphere-Atmosphere Experiment, SCSMEX for South China Sea Monsoon Experiment, KWAJEX 
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Future developments and works 


There is much more work to be done comparing simulated cloud systems over various types of 
land and vegetation environments, ranging from arid to jungle. Recently completed field 
programs (TOGA-COARE, ARM, TRMM LB A, TRMM KWAJEX and NASA CAMEX) 
should provide a good opportunity to orchestrate combined observational and numerical studies 
of convective systems. These large-scale field campaigns can provide some of the desperately 
needed observations for key locations. These observations can guide and correct existing 
microphysical schemes used in the CRMs. 

Recently, physical processes represented in the spectral bin-microphysical scheme has 
been implemented into the two-dimensional version of the GCE model. The formulation of the 
microphysical processes is based on solving stochastic kinetic equations for the size 
distribution functions of water droplets (cloud droplets and raindrops), and six types of ice 
particles: ice crystals (columnar, plate-like and dendrites), snowflakes, graupel and frozen 
drops. Each type is described by a special size distribution function containing 43 categories 
(bins). The bulk density is equal to 0.9 g cnr 3 for ice crystals. Snowflakes, graupel and 
frozen drops are assumed to be spheres and their densities range from 0.01 g cm' 3 to 0.9 g 
cm* 3 . The terminal fall velocities used are those applied by Khain and Sednev (1995), List and 
Schemenauer (1971) and Cotton et al. (1986). Nucleation (activation) processes are based on 
the size distribution function for cloud condensation nuclei (43 size categories). The GCE 
model using the spectral bin-microphysics can be used to study cloud-aerosol interactions and 
nucleation scavenging of aerosols, as well as the impact of different concentrations and size 
distributions of aerosol particles upon cloud formation. These findings will, in turn, be used to 
improve the bulk parameterizations. With the improved GCE model, it is expected to lead to a 
better understanding of the mechanisms that determine the intensity and the formation of 
precipitation for a wide spectrum of atmospheric phenomenon (i.e., clean or dirty environment) 
related to clouds. 

In addition, cloud microphysical processes, heat fluxes from the warm ocean, land and 
radiative transfer processes should interact with each other. How these processes interact 
under different environmental conditions should be a main focus of modeling studies in the 


for The Kwajalein Experiment, CAMEX for Convection and Moisture Experiment, and TEFLUN for 
TExas/FLorida UNderflights Experiment 



future. Also, a major area of needed development involves scale interactions and how cloud 
processes must be included in simulations of mesoscale to global-scale circulation models. 
Specifically, Moncneff and Tao (1999) suggested that improved CRMs can be used to 
addresses the following aspects in the near future: 

1 . Derive physically based parameterizations for numerical weather prediction models and 
climate models; 

2 test s,n g |e - column representations of physical processes (i.e., the processes that trigger 
convection, cloudiness and convective momentum transport); 

J complement large-scale field experiments that would otherwise be subcritical in terms 
of cloud-scale measurements; 

4. add value to data sets in situations where standard soundings are the only measurement 
available; 

5. improve the physical basis of surface (land and ocean)-atmosphere interaction in 
coupled climate models; 

6. help in the design of space-based and earth-based remote sensing and in the 
interpretation of the data sets; and 

7. understand the vortex formation that may be important for initial tropical cyclone 
(hurricane) development. 

Since the real atmosphere is three-dimensional, three-dimensional cloud resolution 
model simulations are also needed to address the above scientific problems. 
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Table Captions 

Table 1 Major foci of cloud resolving model (CRM) development in the past four decades. 

Some (by no means all) key contributors are also listed. (The author apologizes for 
omitting any other major contributors to CRM development.) 

Table 2 Key developments in the cloud resolving model (CRM) approach for studying 
tropical deep convective systems over the past two decades. 

Table 3 Characteristics of the Goddard Cumulus Ensemble Model 

Table 4 Applications of the Goddard Cumulus Ensemble (GCE) Model. The specific topics 
and their respective GCE model characteristics, major results and references are 
shown. 

Table 5 Initial environmental conditions expressed in terms of CAPE, precipitable water and 
Richardson number for the TAMEX, EMEX, TOGA COARE and PRE-STORM 
MCSs. 

Table 6 Values of the ratio R for the GCE model simulations of several convective systems, 
as well as for the 6 h periods corresponding to the initial and mature stages. Their 
respective stratiform rain components are shown in the first column. 

Table 7 The same ratio defined in Table 3.2.2 except for different MCS cases, A, B and C 
of Leary and Houze (1980), cases I and II of Gamache and Houze (1983), the 10-11 
June squall line of Gallus and Johnson (1991), and the COPT squall line of Chong 
and Hauser (1989). Ratios from other CRM studies of convective-stratiform 
interaction are also shown. 

Table 8 Summary of previous cloud-radiation modeling study results. The percentage 
increase or decrease in surface precipitation due to longwave and shortwave effects 
are given along with the mesoscale lifting, if used, for each case. 

Table 9 Ratio of fractional cloud coverage (R = cloud updraft coverage / cloud downdraft 
coverage). Fractional coverage occupied by cloud drafts and active cloud drafts over 



the domain are also shown within the parentheses. This figure is from Tao, 
Simpson and Tao (1987). 

Table 10 The GCE model simulated domain averaged surface rainfall (mm), stratiform amount 
(percentage) and microphysical processes (ratios between evaporation and 
condensation, sublimation and deposition, and deposition and condensation), (a) is 
for TOGA COARE and (b) for GATE. For TOGA COARE, the GCE model results 
are also separated into sub-periods, deep strong convection during December 20-23 
and 24-25 and weaker convection prior to, in between, and after the deep convection 
(December 19-20, 23-24, and 25-26, 1992). Slow-moving (non-squall, September 
2-4), fast-moving (squall, September 4 to 6) and less organized (September 6 to 8) 
periods for the GCE model simulated GATE results are also shown. 
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Figure Captions 

Fig. 1 Estimated surface rainfall intensity over part of horizontal domain at (a) 270 min, 
(b) 285 min, (c) 300 in and (d) 315 min simulated time from a three-dimensional 
GCE model. The contour interval is 10 mm h‘* starting at 1 mm h'l. This type of 
merger is identified as a parallel cells case. 

Fig. 2 Same as Fig. 3.1.1 except for the case of perpendicular cells. Time is (a) 255 min, 
(b) 270 min, (c) 285 min and (d) 300 min. 

Fig. 3 Depiction of trajectory paths calculated from the evolving three-dimensional model 
wind fields for the parallel cell case shown in Fig. 3.1.1. (a) Three-dimensional 
depiction of upward paths as viewed from NNW, (b) as viewed from the overhead. 
They are computed forward from 300 min to 340 min. (c) and (d) are the same as 
(a) and (b). respectively, except they are computed backward from 300 min to 270 
min. The shaded area in (a) and (c) indicates the three-dimensional depiction of 
estimated 20 dBZ iso-surface at 300 min. 

Fig. 4 (a) and (b) are schematic cross-sections through mature squall lines observed in the 

tropics and midlatitudes, respectively. They are adopted from Houze (1977) and 
Smull and Houze (1987). (c) and (d) are the schematic diagrams of major 
conceptual models of tropical and midlatitude squall lines derived from case studies, 
respectively. Areas of light and dark stippling indicate areas of high- and low- 6^ 

air, respectively. They are originally shown in Zipser (1969) and Newton (1963), 
but are adopted from Rotunno et al. (1988). 

Fig. 5 The water budgets for (a) TAMEX, (b) EMEX, (c) PRE-STORM and (d) TOGA 
COARE simulated squall-line MCSs. Italic numbers indicate the amount of 
condensate transfer between various regions and layers while quantities in 
parenthese are the net condensation generated through microphysical processes. 

Fig. 6 Schematic diagram demonstrating the effects by different cloud-radiation 
mechanisms (cloud-top cooling and cloud-base warming - alters the thermal 
stratification of the stratiform cloud layer; differential cooling between clear and 
cloudy regions - enhances dynamic convergence into the cloud system; and the large- 
scale radiative cooling - destabilizes the large-scale environment). 
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Fig. 7 (a) Convective and (b) stratiform heating profiles stored in the heating profile look- 

up table for the CSH algorithm. The profiles were obtained from GCE model 
simulations for cloud systems in various geographic locations [the Pacific warm pool 
region (TOGA1, TOGA2, TOG A3, ERECT - a squall system with erect updrafts, 
MRSH ISL - Marshall Island), the East Atlantic region (GATE), midlatitude USA 
(PRESTORM - PSTM), and Australia (EMEX)]. (c) and (d) are the same as (a) 
and (b) except that these profiles are from Gallus and Johnson (1991, curve GJ), 
Yanai et al. (1973) but partitioned into convective and stratiform components by- 
Johnson (1984, curve Y), Houze (1989, curve H), Houze and Rappaport (1984, 
curve HR), Chong and Hauser (1990, curve CH) and an African squall line 
simulated by Caniaux et al. (1994, curve COPT81). 

Fig. 8 Monthly (February 1998) mean latent heating at (a) 8, (b) 5 and (c) 2 km over the 
global tropics derived from the CSH algorithm. 

Fig. 9 Heating rates for condensation, c, and evaporation, e, in the two-dimensional model 
(dashed line) and the three-dimensional model (solid line). The heating rate estimated 
from the large-scale observations, Ql - Qr, and the total cloud heating rate in the 
three-dimensional model are also included. This figure is from Tao, Simpson and 
Tao (1987). 

Fig. 10 (a) The vertical profiles of the heating rate by condensation of moisture, c, 

evaporation of liquid water drops, e, net vertical flux of the sensible heat, F, the total 
heating rate by clouds and the heating rate estimated from large-scale observations, 
Ql - QR- (b) The vertical profiles of the moistening rate by condensation of 
moisture, evaporation from liquid water drops, net vertical moisture flux, the total 
moistening rate by clouds and the moistening rate estimated from the large-scale 
observations, - (c p /L)Q2. This result was from a 3-D GCE model simulation for a 

GATE convective system (Tao and Soong, 1986). 

Fig. 1 1 Time sequence of the GCE model estimated domain mean surface rainfall rate (mm 
Ir 1 ) for (a) TOGA COARE and (b) GATE. This type of CRM diagnostics and 
graphical presentation have been very popular were first presented in Tao and 
Simpson (1984). 
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Parameters /Processes 

GCE Model 

Vertical Coordinate 
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Explicit Convective Processes 

2 class water & 2 moment 
4 class ice 

Implicit Convective Processes 

Betts and Miller, Kain and Frisch 

Numerical Methods 

Positive Definite Advection for Scalar Variables; 
4-th Order for Dynamic Variables 

Initialization 

Initial Condition with Forcing 
from Observations/Large-Scale Model 

Radiation 

Broad-Band in LW; Solar 
Explicit Cloud-radiation Interaction 

Sub-Grid Diffusion 

TKE 

Two-Way Interactive Nesting 

Radiative-Type (2D model only) 

Surface Energy Budget 

Force-restore Method 
7-Layer Soil Model (PLACE) 
TOGA COARE Flux Module 
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Cloud-Cloud 
Interactions and 
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Ql and Q2 Budgets 


Cloud 

Characteristics 
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Momentum 

Transport 


Ice Processes 


Convective and 
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Interactions 


Cloud Radiation 
Interactions 
& diurnal 
variation of 
* precipitation 


Cloud Chemistry 
Interactions 


Air-Sea 

Interactions 


Precipitation 
Efficiency (PE) 


Model 

Characteristics 


2D/3D 
Warm rain 


2D/3D 

Warm rain and 
Ice Processes 


Major Results 


Cloud downdraft and its associated 
cold outflow play major role in cloud 
merger 


Importance of evaporative cooling in 
Ql budget 

Importance of vertical transport of 
moisture by convection in Q2 budget 


2D/3D 
Warm rain 


2D/3D Identify the role of horizontal 

Smaller Domain pressure gradient force on up-gradient 
in 3D transport of momentum. 


2D/3D 



Active convective updrafts cover 
small area but major contributors in 
mass, Qi and Q2 budgets. 
Excellent agreement with aircraft 
measurements. 


The importance of ice processes for 
stratiform rain formation and its 
associated mass, Qi and Q2 budgets. 


The horizontal transport of 
hydrometeors and water vapor from 
convective towers to stratiform region 
are quantified. 


Longwave cooling can enhance 
precipitation significantly for 
tropical cloud systems, but only 
slightly for midlatitude systems. 

Modulation in relative humidity by 
radiative processes is major reason for 
diurnal variation of precipitation. 


Significant redistribution of trace 
gases by convection. Enhancement of 
O3 production related to deep 
convection in tropics. 


TOGA COARE flux algorithm 
performs well compared with 
observation, better than other flux 
algorithms. 

Surface fluxes are important for 
precipitation processes and maintain 
CAPE and boundary layer structure. 


Examined different definitions of PE. 

Identify several important 
atmospheric parameters for better PE. 


2D (short and 
long term 
integration) 


2D/3D 
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Tao and Simpson 
(1984, 1989a) 
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( 2000 ) 
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Tzioetal (1987) 


Soong and Tao (1984), 
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al (1991), Tao et al . 
(1993a), Ferrier et al 
(1995) 


Tao et al. (1993a), Sui 
et al (1994), Tao 
(1995), Lang et al. 
( 2000 ) 
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Tao et al. (1996), Sui 
etal (1998) 


Thompson et al (1997 
- a review) 


Wang et al (1996, 

2000) 


Ferrier et al. (1996) 


Table 4 


































Land Processes 

2D/3D 

Importance of mesoscale circulation 
induced by soil gradient on 
precipitation. Identify the 
atmospheric parameters for 
triggering convection. 

Lynn et al (1998, 
2000a, b). Baker et al 
(2000) 

Idealized Climate 
Variations in 
Tropics 

2D 

Examined several important 
hypotheses associated with climate 
variation and climate warming. 

Identified physical processes that 
cause two different statistical 
equilibrium states (warm/humid and 
cold/dry) in idealized climates. 

Sui et al (1994), Lau 
et al (1994, 1995) 

Tao et al (1999, 
2000), Shie et al. 
(2000) 

TRMM Rainfall 
Retrieval 

3D 

Improved the performance of TRMM 
rainfall retrieval algorithms by 
providing realistic cloud profiles. 

Simpson et al (1996 - 
a review) 

Latent Heating 
Profile Retrieval 

2D 

Developed algorithms for retrieving 
four dimensional vertical structure of 
latent heating profiles over global 
tropics. 

Tao etal. (1990, 
1993b, 2000a, b) 
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